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ABSTRACT 


Coming to periscope depth 1s one of the most intensive of the routine submarine 
operations. Errors in Fire Control and Sonar System information serve to produce 
uncertain contact solutions that complicate the decision of selecting a safe course. The 
model developed in this thesis simulates a specified number of trials on each possible 
course, with the measure of effectiveness for each course being the probability of the 
course being acceptable with respect to specified minimum range criteria. The model 


outputs a geographic display and a graph of the measures of effectiveness versus course. 
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EXECUTIVE SUMMARY 


Coming to periscope depth is one of the most intensive of the routine submarine 
operations. Errors in Fire Control and Sonar System information serve to produce uncertain 
contact solutions that complicate the selection of a safe course. The model developed in this thesis 
simulates a specified number of trials on each possible course, with the measure of effectiveness 
for each course being the probability of the course being acceptable with respect to specified 
minimum range criteria. The model outputs a geographic display and a graph of the measures of 
effectiveness versus course. 


Specifically: 


1. The model determines all safe courses with respect to either of two user-specified 
range criteria using a simulation model. The model is in a readily available 
programming language, Excel, with minimal yet thorough manual inputs from available 
data sources on board. This course scoring includes constraints pertaining to the 
current and time-advanced contact situation, the errors associated with each term 
defining the contact situation, and the directions of seas. Data inputs are limited to 


realistic information available for a timely ascent. 


2. The model formulates a standard tactical graphical presentation depicting the scoring 


of all available courses. This visual aid displays: 
a) range and bearing to each contact, 
b) range and bearing probability ellipse for each contact, 
c) the 30 minute dead reckon true track for each contact, 
d) the scoring of each course for the given scenario, 


e) atable format summary of the input information. 


1X 


3. The model provides a visual aid displaying the measures of effectiveness for each 
course. The visual aid also displays the preferred course sector for the given direction 


of seas. 


I. INTRODUCTION 


One of the most perilous of the routine submarine shiphandling operations is 
proceeding to periscope depth. This operation 1s a prelude to several evolutions, most 
notably surfacing. A course must be decided upon that will ensure safe completion of the 
ascent. A safe course is any course with an acceptably low collision risk, essentially zero, 
with any surface contact. A safe course must keep all contacts in excess of some specified 
minimum range, and it must be selected expeditiously. Three main factors contribute to 
the level of danger: 

1. inaccuracies of Fire Control solutions and Sonar information, 

2. own ship’s presence and intentions are unknown to surface contacts, 

3. environmental conditions. 

While several sources of information are available, both equipment and operator 
induced inaccuracies must be considered. Collection of further information with own 
ship’s maneuvers for solution refinement may be limited due to a need for a timely ascent.’ 
Perhaps most dangerous, environmental conditions may contribute to partial or complete 
masking of contacts and limit the range of feasible courses. In areas of very high contact 
density, a submarine may have as many as ten significant contacts to account for during 
the ascent. To relate the difficult nature of this maneuver, one could imagine landing a 
plane with no windows in an active parking lot full of vehicles unaware of the plane's 
descent. Additionally, each vehicle present will have the nght of way with respect to the 
approaching aircraft. 

The objective of this thesis 1s to develop a tactical decision aid to assist the Officer 
of the Deck in determining the acceptable courses upon which to proceed to periscope 
depth. An acceptable course will place own ship clear of any contacts or other 
navigational hazards. Additionally for enhanced ship control, courses in a + 30° sector on 


either side of the approaching seas are preferred. 


1. Emergency ascents will not be covered. 


The next chapter provides background information relevant to the problem and the 
model development. Chapter III describes the methodology for the problem formulation 
including any pertinent assumptions. Chapter IV addresses the development of the model. 
Chapter V provides an overview of the analysis performed on the output data with respect 


to variations of the input information. The final chapter states the conclusions. 


Hi. ASCENT TO PERISCOPE DEPTH 


A. BACKGROUND 


A submarine comes to periscope depth for a variety of reasons in addition to 
performing a visual search. Operations at periscope depth are necessary for receipt of 
communications, navigational fixes, ventilation of ship, and as a prelude to surfacing. 

In all but an emergency ascent or surfacing, some preliminary maneuvering is 
required to search for yet ungained contacts concealed by the baffles or the thermal layer. 
The baffles are the bearings along which own ship’s sonar equipment has reduced 
capabilities due to the location of the sonar sensor on own ship. For example, a forward 
mounted sonar array will not detect contacts in a given sector aft of own ship. The baffles 
are much like a driver’s blind spot. Unlike checking the mirrors when changing lanes, a 
quick glance and some additional speed are not sufficient. The effects of a thermal layer 
will be covered in Section 2.a below. 

Own ship’s maneuvers are used to resolve the essential parameters which define 
the motion of the contact(s). The process by which these parameters are determined is 
known as Target Motion Analysis, TMA. By use of judicious maneuvers and the 
conservative assumption of a closing contact, the essential information upon which to base 
a periscope depth course decision can be obtained in a timely manner. As critical as this 
periscope depth evolution is, only a few basic parameters about each contact are required. 
While several techniques are available for refining the contact data to determine the 
greater details of its motion, these greater details are not required to determine an 
acceptable course upon which to proceed to periscope depth. Consistent with a 
submarine’s need to remain undetected, only passive TMA techniques are employed. 

Ship’s speed is limited in the ascent to periscope depth due to the hydrodynamic 
force exerted on the extended periscope mast and fairing. This speed limitation also 


somewhat restricts own ship’s maneuverability due to the reduced steerageway. This 


effect is due to the change in the amount of hydrodynamic force exerted on the same 
rudder area at lower ship’s speeds. The reduced maneuverability serves to intensify the 
evolution. It should be noted that preascent maneuvering does not have this same speed 
limitation. 

Now that the basic problem has been stated, the factors affecting the periscope 


depth course will be addressed. 


B. FACTORS 


1. Contact Situation 


a. Contact Motion Parameters 


A contact is any detected underway vessel. When considering possible 
periscope depth courses, each contact’s motion must be considered. The essential 
parameters are: 

1. relative position, determined by range and bearing from own ship, 
2. true motion, determined by course and speed of each contact. 
Determination of these parameters is possible through Bearing-Only TMA 


or from the Fire Control and Sonar Systems. 
b. Associated Errors 


The major problem with using the information from these latter sources is 
the errors associated with the generation of the information. These errors are both 
intrinsic in the equipment, be it from physical design or analytical methods used, and 
operator-imposed through human error, interpretation, or level of experience. The 
method by which these errors will be accounted for will be through simulation modeling 
focused on determining the effects of such errors and a probabilistic scoring method to 


display the viability of any available courses in the presence of these uncertain errors. 


2. Environmental Conditions 


In the broad spectrum of environmental conditions that affect submarine 
operations, only two significantly affect the periscope depth course chosen - one directly 


and one indirectly. 
a. Thermal Layer 


The presence or absence of a thermal layer can impact the contact scenario 
by concealing ungained contacts until the final ascent. Due to refraction, sound waves 
traveling through distinct transitions in density of the given medium, in this case seawater, 
are deflected or bent away from the regions of lower density, or warmer water. Therefore 
a warmer layer of seawater below the surface will cause surface noise to bend back 
towards the surface and sounds generated below the layer to be deflected towards the 
ocean bottom. As a result, sensors below the layer will not detect surface noise. The 
stronger the thermal layer, the greater the potential for gaining new contacts on ascent. 
The error terms associated with such contacts are greater due to the lack of solution 


refinement. 
b. State and Direction of Seas 


The direction of seas is the direction from which the seas are approaching. 
This direction is determined by acoustic trace patterns on certain sonar displays. Similar 
aural techniques can be used to determine sea state, an important factor in deciding how 
much relative weighting should be given to the direction of seas. The higher the sea state 
the more prominently shiphandling will be affected. It is preferred to come to periscope 
depth with the direction of seas within 30 degrees of the bow. Courses in this sector aid in 
ship control by maximizing the relative speed between own ship and the local current. 
This not only assists in maximizing steerageway, but also limits the two undesirable effects 


described in the next two paragraphs. 


If the chosen course is perpendicular or nearly perpendicular to the 
direction of seas, excessive roll may be experienced, potentially limiting ship control. 

If the chosen course is along or nearly along the direction of seas, ship’s 
control is affected due to reduced steerageway. More importantly, with the seas coming 
up the stern, the relative speed of the seas with respect to own ship’s speed is reduced. 
Because of this, own ship is more susceptible to the wave motion. As the waves travel 
along the longitudinal axis of the boat, depth control becomes difficult due to the 
“porpoising” effect caused by the resulting undulating pitch motion. 

These effects are amplified in higher sea states and are next to negligible in 


calm seas. 
3. Geographic or Navigational Constraints 


Geographic or navigational constraints will not be considered in this model, as they 
are easily accounted for in the course decision. The option does exist to allow entry of a 
navigational hazard as a contact with no speed yet still possessing the error terms 


associated with its range and bearing. 


I. METHODOLOGY 


The methodology will be addressed from the standpoint of a single simulation trial. 
The model simulates over all possible own ship’s courses to determine the acceptability of 
each course with respect to each specified range criterion encountered over a given time 


on course. 
A. INPUT DATA 


The simulation requires manual entry of all parameters. This entry is via dialog 


boxes. Specifics of the entry methods will be addressed in Chapter IV. 
1. Required Input 


The required parameters are categorized according to their use within the model. 
The parameters are: 
IL. Scenario Parameters 
a) NumberOfContacts, Number of Contacts to be entered, 
maximum of five 
b) Rhccminrn, Acceptable Minimum Range to any contact 
during time on course, in yards 
c) Rhsarery, Safety Range to any contact during time on course, 
in yards, always less than Rhyccuawry 
d) TimeOnCourse, Time on Course, in minutes 
e) DMho, Own Ships Speed, in knots 
w Simulation Parameters 


a) NumberOfTrials, Number Of Trials per course to be run 


b) StepSize, analyze every n™ course, integer value from one to 
seven 
3. Contact Motion Parameters (for each contact) 


a) Ct, Target Course, in degrees 

b) DMht, Target Speed, in knots 

Cc) Rh, Range, in yards 

d) By, True Bearing to Contact, in degrees 
e) CtSigma, Course Error, in degrees 

f) DMhtSigma, Speed Error, in knots 

g) RhSigma, Range Error, in yards 

h) BySigma, Bearing Error, in degrees 


B. TARGET MOTION ANALYSIS 


All Target Motion Analysis is performed in the relative frame of reference. This 
accommodates determination of the range and time of closest point of approach as well as 
direct use of error terms in the simulation. All TMA is performed with own ship and all 


contacts on constant course and speed. 
1. Translation between Coordinate Systems 


Translation from the polar coordinate parameters of range and bearing is 
performed to allow use of Cartesian coordinates both in the determination of time of 
closest point of approach, fcp, and the minimum range encountered, R/ayy. All 
computations are performed using the standard mathematical axes versus the standard 
tactical axes. Doing so simplifies the required computations within the simulation. The 
translation to the standard tactical axes is performed only as required to generate the 
geographic display. 

The translation from polar to Cartesian coordinates 1s: 

X pos = Mh: cos{ By) 
ee ie sin( By) 
Similarly own ship’s and target’s course and speed are translated to Cartesian 


coordinates. Co is own ship’s course. 


X, = DMho-cos(Co) 
¥, = DMho-sin(Co) 
X,, = DMht -cos(Ct) 
Y¥. = DMht -sin(Ct) 

These conversions to Cartesian coordinates also allow the direct computation of 
points necessary to produce the dead reckoned traces present in the Geographic Display 
Graph. 

The input parameters of DMho, Ct, and DMht along with Co as generated by the 


simulation are used to determine the Cartesian components of relative motion, xDMhr and 


yDMhr. 





yds 
DMhr =| X,- X_ }kt (3375 
rai ( . o) : kts- min 


yds 
DMhr =\Y. — Y.. tkt ( a. on 
y { (¥, 0) 8 2 aes 


2. Incorporating Error Terms into Target Motion Analysis Parameters 


a. Error Term Distributions 


Each of the four TMA parameters, Rh, By, Ct, and DMht;, have an 
associated input error term, RhSigma, BySigma, CtSigma, and DMhtSigma respectively. 
Each error term 1s used as twice the value of a standard deviation in that parameter. The 
mean of each error distribution is assumed to be zero. The distributions of the error 
terms used by the simulation are independent Normal ( 4 = 0, o =stated error/2). 

It 1s critical that the operator realize the mathematical transition of the 
stated error to the standard deviation of the error term distribution. The probability that a 


normal error falls within a stated + interval varies as shown in Figure 1. 


Probability in + Interval 





Figure 1. 

These probability values were found using the standard normal tables and 

the cumulative distribution function as follows: 
p=P(X 2ko)- P(X <—-ko), for k =1,2,3. 

These values demonstrate the judgmental subjectivity with manual 
assignment of the error values. One operator may state a + interval gauged from the 
perspective that the interval will contain the contact 90% of the time whereas another 
Operator may gauge that percentage differently. For the purposes of this thesis error 
values are to be expressed as two standard deviations. This 1s to say, the specified contact 
parameter will be within the error bounds with probability 0.9544. The input dialog 
includes a note above the error term entry boxes to this effect. 

The normal distribution generates a probability density according to: 

eee 


f(x) = —— 


The cumulative distribution 1s therefore: 








As no closed form solution exists for the normal cumulative distribution, 
solving for x directly is not possible. Instead, a value of the cumulative probability, /(x), is 
supplied and a numerical method is applied to attain a value of x. An iterative normal 
inversion function, NormInv(p = probability, 4 = mean, o = standard deviation), 1s 
available in Excel, but it is expensive in regards to computational time. As an alternative 


to calling NormInv(y, 4, o) several times within each tral, the standard normal was 


10 


discretized into 21 equally sized bins with the values of the standard normal curve being 
assigned to the elements of an array. In doing so, each bin has a value equal to the 
number of standard deviations from the mean for the equally sized steps of cumulative 
probability from 0 to 1 by increments of 0.05. The values of 0 and 1 are avoided to 
prevent faulty returns from the normal inversion function call. The values of 0.001 and 
0.999 are used instead. 

The curve, (x), used to generate the bin values is plotted below in Figure 
2. The vertical axis is the uniform random probability. The horizontal axis is the number 


of standard deviations from the mean. 


random probability 





Standard Deviations from the Mean 


Figure 2. 

A uniform random number is generated, multiplied by the number of bins, 
and truncated. This truncated value is the index of the array element to be used in 
determining the error term. The array element is multiplied by the standard deviation, o. 

p = random number 
index = int(p * bins) 
error = array(index)*o 

The resulting value is an approximately normal error term with a mean of 
zero and a standard deviation of the stated error, 20 divided by 2. 

The resulting error terms for Rh, By, Ct, and DMht are RhError, ByError, 
CtError, and DMhtError respectively. These values are regenerated with each trial and 


incorporated as follows: 


11 


trialX , = (Rh+ RhError)- cos( By + ByError) 
trialY, =(Rh+ RhError)-sin( By + ByError) 
trialX , = (DMht + DMhtError) -cos(Ct + CtError) 
trialY, = (DMht + DMhtError) - sin(Ct + CtError) 
The trial values of the parameters are sent to the subsequent function calls 


for determination of the closest point of approach. 
b. The Independence Assumption 


To justify the use of independence between the associated error terms, a 
dependent case example is presented as a sensitivity test. In passive TMA, contact bearing 
rate determines the contact speed across the line of sight rather accurately, [Ref. 2.]. This 
value is key to determining Ct and Dmht. Any DMhtError induces a corresponding 
CtError in the dependent case as DMht is used in attaining Ct or vice versa. As passive 
TMA may be the sole technique used by a submerged submarine, this case of (Ct,DMht) 
dependency is of interest. The dependent version of the model applies dependence 
between Ct and DMht as follows: 


Let xDMht be contact speed across the line of sight, assumed known 

Let depDMAht be the dependant DMht 

xDMht = DMht - sin( By — Ct) 
sin(trialBy —trialC t) 

It should be noted that this case of dependence is one of many possible 
dependent relationships. This particular case is chosen for its straightforward nature. The 
results of the dependent version of the model will be addressed 1n Chapter V. 

The input error terms for the simulation model are an instantaneous 


evaluation of the contact’s solution accuracy as of the input time. The model makes no 


attempt to modify the magnitude of any error terms over the specified 7imeOnCourse 


parameter. This is consistent in that the model does not account for any solution 


refinement or lack thereof during the time on course considered. 
3. Minimum Range Determination 


The time of closest point of approach is the time at which the relative motion 
vector is perpendicular to own ships position. To determine the time at which CPA 


occurs, the range function is minimized. 
RA(t) = V(X p05 + XDMhr - 1)” + (Ypos + yDMhr - t)’ 
To find the time at which range is a minimum, the first derivative with respect to 


time is taken. 


= (Rh(0) 7 “ Cag PDE Ce + yDMhr - 1)" 


d (Rh ) 2| xDMhr (X05 + xDMhr -t) + yDMhr - (Yp95 + YDMhr - t)| 
_— LQ) = ——————S——— 
dt 24/(X pos + XxDMhr - 1)? + (Vong + yDMhr - 1)” 

By setting the first derivative of the range function to zero and solving for time. 


xDMhir «(X poe + xDMhr - t) + yDMhr (Yn, + yDMhr -t) = 0 
xDMhr - X po, + xDMhr* -t + yDMhr -Y,5, + yDMhr* -t =0 


~ xDMhr - X ppg — DM - Yves 


"PA xDMhr? + yDMhr? 

The geometry of the problem dictates this value to be a minmum. Next, fcp, 1s 
used to determine the minimum range experienced to the given contact over the 
TimeOnCourse. 

If tcp, iS negative, or equivalently, before the time interval specified by 


TimeOnCourse, the contact is now opening in range; therefore the minimum squared 


range 1s: 
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sla = (RA(0)y 
If tcp, is within the time interval as specified by TimeOnCourse, then the minimum 


squared range is: 


iy = (Rh (1 CPA y 
If tcp, is greater than 7JimeOnCourse, then the contact closes in range until the end 
of the time interval, so 
Rh?,, = (Rh(TimeOnCourse)) 
This illustrates an important consideration in using this model in that it does not 
account for contact ranges at times beyond that specified by the Time on Course 


parameter. 
C. MEASURES OF EFFECTIVENESS 


The measures of effectiveness, MOLEyccumrn and MOE surgery , are the percentage 
of trials for which all contacts remain at ranges in excess of the acceptable minimum range 
or safety range on any given course. Equivalently stated, each MOE is the probability that 
the course is acceptable as defined by not violating the applicable range constraint. The 
squared range is again used. 

The first step in determining the MQOEs is to take the minimum of the vector of 
minimum squared ranges (the use of squared ranges avoids the square root function) to all 
contacts as follows: 

Let Rh? be the minimum squared range to contact i, and 
a min} Rh. Fee ee 

If z is greater than the acceptable minimum range squared, then n, a counter 
variable, is incremented. If the value of the minimum is greater than the safety range 
squared, then 4, another counter variable, is also incremented. After all simulation trials 


have been run, the MOEs are calculated as follows: 
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Let n f be final value of n 
Let k, be final value of k 





a 
MOE sccuavna = trials 
ky 
MOE sarery = trials 


MOE accminen 18 then the fraction of trials for which all contacts remain in excess 
of Rhaccumnrn, and MOE surery is the fraction of trials for which all contacts remain in 
excess Of Rhsirery. The value of each MOE is then stored in an internal array and 


displayed on the worksheet “Graphs”. 


D. OUTPUT GRAPHS 


Two distinct output graphs are produced to provide not only a display of the 
MOEs of each course, but also to visually verify the results against a geographic display of 


the specified contact scenario. 
1. Measures of Effectiveness Graph 


Each of the MOEs is plotted as a scatter graph versus the entire range of courses, 
0 to 359. The MOEsccvien is plotted in front of the MOEssrery. This is possible as the 
value of MOE sccumven is always less than the value of MOEsurerr, given correct inputs. 

The use of scatter plotting assists in showing the acceptable course sector(s) 
regardless of the specified StepSize. The graph also displays the courses allowed by the 
safety range constraint in the event that few or no courses meet the acceptable minimum 
range constraint. Additionally, spikes displaying the + 30° preferred course sector with 


respect to the direction of seas are included. See Appendix A. Output #1, 2, and 3. 
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2. Geographic Display Graph 


The geographic display was created to provide for internal verification of the 
model. The graphs provide an overview of the contact’s position, motion, and the 
probability ellipse for RA and By errors. First an ellipse with axes of Rh error and By error 
is created in Cartesian coordinates. The ellipse is then rotated by an angle equal to the 
contact’s bearing to align the ellipse such that the RA error axis is along the radius of a line 
extending from the center out in the direction of the contact’s bearing. The ellipse is then 
translated out to the contact’s Rh and By through a conversion to Cartesian coordinates. 

In addition to the (Rh, By) probability ellipse, the contact’s course, Ct and speed, 
DMhi are reflected by dead reckoning each contact in two minute intervals out to 30 
minutes. This permits a visual indication of the contact’s motion. It does not include any 
visual representation of the Ct and DMAhrt errors. 

The Geographic Display Graph provides not only for viewing the contact scenario 
on a standard tactical display but also shows the scoring of all courses with respect to the 
acceptable minimum range criteria as a value proportioned to the maximum initial range to 


any contact. This method shows points created by the following: 
Score = max{Rh, (0),...,RDyumseropconane (0)}* MOE goenanet 
Using this relationship, scores at the maximum range of any contact have 
MOE accuirn = 1, at 50% of maximum range, MOEyccumvrn = 0.5, etc. Similar scoring 


points for safety range course scoring were excluded so as to not visually overload the 


display. See Appendix A, Examples #1, 2, and 3. 
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IV. SEMULATION MODEL 


A. DESIGN 


This section will address the various design elements of the model. The model 
consists of approximately 410 lines of code and one worksheet to support the graphing 


features. 
1. Programming Language 


The model is coded in Visual Basic for Applications, Windows Version 4.0, within 
Excel Version 7.0 of the Microsoft Office 95 Suite. This programming language was 
chosen due to extensive versatility with respect to mathematical, statistical, and graphical 
display capabilities. 

The choice was also in consideration of availability of the underlying program in 
the event the model is deemed useful as either a tactical decision aid or a training aid. Use 
of this readily available program allows any user access to the model with no additional 
programs or system capabilities beyond those of an average PC, once a copy of the 
program has been provided. Appendices B, D, and E contain the coding for the input, 
simulation, and basic output. Recreation of the graphs is left to the user or will be 


provided upon request. 
2. Modularity 


The source code in separated into three modules, “Data Entry”, “Simulation”, and 


“Main”. 
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a. Data Entry Module 


The “Data Entry” module includes the procedures required to create the 
graphical user interfaces, GUIs, necessary to allow the user to enter the required 
parameters. This module also contains all the control procedures for the GUIs. The 
initialization procedure is also contained within this module. Lastly, the procedure that 
sends the input parameters to the “Graphs” worksheet is located in this module. The 


“Data Entry” module code 1s presented in Appendix B. 
b. Simulation Module 


The “Simulation” module includes the simulation procedure and all the 
functions called by it. Also included in this module is the procedure to send the MOE 
matrix to the “Graphs” worksheet. The “Simulation” module code is presented in 


Appendix D. 
c. Main Module 


The “Main” module contains only the main procedure. The main 
procedure calls the modules that control each major portion of the model; data entry, 
graph inputs, simulation, and graph output. The “Main” module code is presented in 


Appendix E. 
B. IMPLEMENTATION 


1. Control Flow Path 


Control begins with the implementation of the “Main” procedure. “Main” calls 
“PrepareDialogs” which in turn calls “InitializeContacts” and also initializes the edit boxes 


of both input dialog boxes. The procedure then displays each dialog in turn retrieving the 
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data input by the user. The input data is assigned to the applicable variable not by 
“PrepareDialogs” but instead by the individual control’s “Change” procedure. 

After this the control returns to “Main”. Next “Main” calls “SendInputToGraphs’”, 
which sends the input to the applicable cells within the “Graphs” worksheet for use in 
generation of the graphs. Control now returns to “Main”. 

Next “Main” calls “Simulate” by assignment to the variable “Data”. The inner 
workings of the “Simulate” function will be covered in section 3 below. After the 
simulation has ended control returns to “Main”. 

Finally “Main” calls ““SendDataToGraphs”, which sends the MOE matrix to the 
“Graphs” worksheet for use in displaying the applicable course scores on each output 


display. Control then returns to “Main” and the program is terminated. 
2. Input Dialogs 


The “PrepareDialogs” procedure displays each dialog as called. The first dialog is 
contained in the “GetInitialData” dialogsheet. This dialog collects the scenario and 
simulation parameters. See Appendix C, Figure 1. 

The second dialog is contained in the “GetContactData” dialogsheet. It retrieves 
the contact motion and error parameters for each of the specified number of contacts. See 
Appendix C, Figure 2. 

It is important to note that the edit boxes require strictly numeric entries. Any 
transition of the cursor between edit boxes must be performed by exclusive use of the Tab 
key. Movement within the contact number dropdown list of the second dialog is more 


flexible. 
3. Simulation 


The simulation begins by initializing the decision variable, Co. Next, the specified 
number of trials is performed. Within each tral, each contact motion parameter of the 


current contact 1s modified by an error. With these new contact motion parameters the 
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“GetTimeToCPA” function ts called. The returned value of fcp, is used in the function call 
to the “GetMinSquaredRange” function. The value returned from this call is the squared 
value of the minimum range achieved to the current contact. This value is assigned to a 
dynamic array of length commensurate to the specified number of contacts. This process 
is repeated within each tnal for the specified number of contacts. 

After all contact’s minimum squared ranges have been determined, the absolute 
minimum squared range for the current trial is determined by taking the minimum of the 
minimum squared range array. The value of the absolute minimum squared range is now 
compared to the specified values of Rhuccumwry squared and Rhsirgry squared. If the 
absolute minimum squared range is greater than or equal to the applicable range squared, 
the appropriate counter variable is incremented. This entire process is repeated for the 
specified number of trials. 

After all the trials have been executed for the current value of Co, each counter 
variable is divided by the specified number of trials then the value is assigned to the MOE 
matrix being positioned by column according to the applicable range criteria and by row 
according to the current value of Co. 

The Co value is then incremented and the entire process repeated for all remaining 


courses. The final step is to assign the OE matrix to the value of the function variable. 
4. Output Graphs 


The form and formatting of each output graph exists on the “Graph” worksheet. 
Each run of the model updates the graphs to depict the most recent execution of the 
model. This method reduces execution time by not recreating the graphing forms or 
formatting with each execution of the model. It is possible to view the visual updating of 


the graphs once the simulation portion of the model has been executed. 
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C. TRANSPORTABILITY 


1. Add-In Conversion 


The model could be transitioned into an Add-In format for internal inclusion in any 
version of Excel 7.0 and beyond. Add-Ins also afford greater security by allowing 


password protection of the program. 
2. Platform Performance 


The model is capable of reasonable execution times on an average PC. As a 
reference, the model can execute a five contact 100 trial simulation in approximately 10 
minutes on a 120Mhz machine. This is clearly not within the range required to make an 
expeditious course decision on a rapid ascent, though is usable in other scenarios. This 
time will decrease, of course, on faster PC’s. 

It is clear through step-oriented call tracing that actual simulation takes about 85% 
of the total execution time. The remaining execution time is used in writing values to the 
“Graph” worksheet. Although graphs can be created directly from the source code, the 
execution time is increased in the creation of all the graph formatting. The process of 
writing values to worksheets seems to take an undue amount of time. The data transfer to 


the worksheet 1s an area to explore in attempting model improvement. 
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V. ANALYSIS 


The simulation will be analyzed with the use of three examples. These examples 
are not very realistic in that the contacts are arranged symmetrically and it most cases are 
assigned similar motion parameters. This was done to assist in displaying the model’s 
characteristics and in no way represents any limitation of the model. The outputs are 
provided in Appendix A. 

The simulation model is verified through use of reference problems from the 
Maneuvering Board Manual, Pub 217, [Ref. 1.], and through vamnations of all input 
parameters to ensure compliance with anticipated effects and current tactical guidance. All 


verification of the model is internal. No external agencies have assisted in the verification. 


A. VIABILITY OF MEASURES OF EFFECTIVENESS 


In Examples 1,2, and 3, compliance 1s evident between the range of acceptable 
courses shown by the MOE graph and the intuitive course decision based on choosing a 
course which has the contacts on the left, drawing left, and contacts on the night, drawing 
right. Additional courses are shown to be acceptable by the MOE graph as the simulation 
course decision does not prevent use of courses that permit any contact to cross the bow 
of own ship. Use of courses that do not allow any contact to cross the bow of own ship 
are always preferred. These courses provide a greater safety margin to a worsening 
situation caused by contact maneuvers towards own ship. The additional courses shown 
by the ZOE graph are useful in scenarios which have limited courses based on the bearing 


rate method alone. 
B. ERROR VALUE INFLUENCES 


The error values used for each scenario are displayed in the summary table of each 


example. The geographic display shows the (Rh,By) probability ellipse based on these 
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values. While no visual representation of Ct or DMAt errors 1s available, the influence of 
all the error terms is evident in both the scoring ring of the geographic display and on the 
MOE graph. Example #1 in Appendix A most clearly displays the effects of the magnitude 
of the error terms. 

Contact 1, bearing 060°, has the lowest magnitude of error. The influence of the 
error terms is evident by inspection of the edges of the scoring ring in the geographic 
display, and in the edges of the scatter graphs of each MOE in the vicinity of each 
contact. The edges in both displays are steeper and have less variation along their slope 
nearer to Contact 1. The edges in the vicinity of Contacts 2 and 3 are much less steep and 
show greater variation in the scoring ring and the scatter graph slopes. Also the peaks on 
either side of lower error term contacts are higher and the valleys are lower. This verifies 
the correct inclusion of the error terms in that the greater the magnitude of the error 
terms, the less distinct the MOE graph results are for the given situation, as is evident due 
to the symmetry of the contact situation. Example #3 also clearly displays the effects of 
the magnitude of the stated error terms. 

Example #2 is highly asymmetrical. The asymmetry is in support of the discussion 
in sections C and E below. The range of acceptable courses from 95° to 200° concurs 
with the author’s intuitrve course selection made through inspection of the contact 
scenario alone. 

Example #3 displays the justification of the independence assumption of C7? and 
DMht. Contacts 1 and 2, shaded in the summary table, are formulated with the dependent 
version of the model coding, see Section B.2.b of Chapter III. Contacts 3 and 4 are 
formulated using the standard independent version. The RhError and ByError for all 
contacts are set to zero to isolate the effects of the Ct and DMht dependency. The 
asymmetry in the MOEs and scoring ring about the North-South axis displays the effects 
of the independence assumption. In the MOE graph, the dependent results have flatter 
peaks and valleys though steeper edges. The independent case is conservative with 


respect to the dependent in that a more narrow range of courses 1s acceptable for contacts 
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with the same error terms. The use of independent error term distributions also simplifies 
computations within the inner-most loop of the simulation model, thereby reducing 


execution time. 
C. SAFETY RANGE VS. ACCEPTABLE MINIMUM RANGE 


The inclusion of the safety range MOE is valuable in that it allows visual 
representation of any additional acceptable courses with respect to the lower of the range 
criteria. This element of the model is also very useful in showing courses acceptable with 
respect to the safety range criterion when no courses meet the acceptable minimum range 


criterion. Example #2 demonstrates an instance of this case about the course of 265°. 


D. CONTROL OF ACCURACY AND EXECUTION TIME 


1. Accuracy 


The values of the simulation parameters dictate the accuracy of the output. A 
higher NumberOflrials per course results in greater accuracy in the MOE’s. A lower 
value of StepSize results in more thorough information by analyzing a greater portion of all 


available courses. 
2. Execution Time 


Execution time is roughly proportional to the values of NumberOfTrials and 
NumberOfContacts, and inversely proportional to StepSize. The greater the 
NumberOfTrials and/or the greater the NumberOfContacts, the greater the execution 
times. The greater the StepSize, the lower the execution times. The user can select values 
for NumberOffrials and StepSize. The NumberOfContacts is limited to a maximum of 5 
for execution time concerns. The user must consider the tradeoff between the exclusion 
of any contact and the lack of consideration, by the model, of that contact over the entire 


TimeOnCourse. If contacts are to be screened for entry, the contact must be deemed not 
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a concern for the entire 7imeOnCourse for any choice of own ships course. Reference 


values of execution times on a 120MHz PC are as follows: 


Ca ara STEPSIZE | NUMBEROFCONTACTS EXECU a TIME 





Example #2 demonstrates the difference in the WOE graph display for a StepSize 


other than 1. 
E. DIRECTION OF SEAS 


The model does not make any attempt to mathematically incorporate the direction 
of seas into the simulation. The relative weighting of direction of seas with respect to 
course selection is dependent on several factors, most notably the sea state. Due to the 
very subjective nature of the incorporation of the direction of seas into the course 
decision, the model only displays the preferred course sector with respect to the direction 
of seas. Example #2 displays a case in which several courses within the preferred sector 
are acceptable with respect to the Rhg4rery criteria, while none are acceptable with respect 


to the RAyccywirn criteria. 
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VI. CONCLUSIONS 


A. MODEL PURPOSE 


The model provides useful tactical information given standard and readily available 
inputs. Reference problems and multiple challenging tnal scenarios demonstrate 
compliance with the pnmary course selection criterion. Furthermore, the model clearly 
displays courses acceptable with respect to the specified range criteria in a single display, a 


valuable resource not currently available. 
B. THESIS OBJECTIVES 


The simulation model provides a viable base from which to continue development. 
The proposed tactical decision aid is useful in many scenarios in addition to the periscope 
depth scenano. The need for further development lies mainly in reducing execution times 
so as to allow use of up to 1000 tnals per course thereby providing greater accuracy and 


reliability, very necessary qualities in data used towards this level of decision. 
C. TACTICAL USE 


The pnmary inhibitor in immediate tactical use is external validation. Accreditation 
of the model would need to follow. The attractive feature with respect to the 
development is that once these processes are complete, the transition to onboard use is 


expeditious as the model is developed for use on a PC. 


D. TRAINING USE 


The model has potential for use as a training aid, where the problem of execution 


times is less restrictive. The development of a random contact generator would greatly 


2g 


increase the value as a training aid by allowing the user to run more scenarios in the same 


amount of time. 


While the model was developed primarily to focus on the periscope depth issue, it 


is potentially useful in any surfaced or submerged course selection scenario. 
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APPENDIX A. COMPOSITE OUTPUT DISPLAYS 
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APPENDIX B. VISUAL BASIC SOURCE CODE FOR DATA ENTRY AND 
CONTROLS 


Type Contact 
Ct As Single 
CtSigma As Single 
DMht As Single 
DMhtSigma As Single 
Rh As Single 
RhSigma As Single 
By As Single 
BySigma As Single 
End Type 
Type Parameter 
NumberOfContacts As Single 
AcceptableMinRange As Single 
SafetyRange As Single 
TimeOnCourse As Single 
DMho As Single 
DirectionOfSeas As Single 
NumberOfTnials As Single 
StepSize As Single 
Bins(0 To 20) As Single 
End Type 
Type CCR 
XO As Single 
Yo As Single 
End Type 


Public Contacts(1 To 5) As Contact 


Public Parameters As Parameter 


Public CCRider As CCR 
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Sub InitializeContactsQ 
Dim Count As Single 
For Count = 1 To 5 

With Contacts(Count) 
(eu 0 
.DMht = 0 
Wn et) 
By =0 
.CtSigma = 0 
.DMhtSigma = 0 
.RhSigma = 0 
.BySigma = 0 

End With 

Next Count 
End Sub 


Sub PrepareDialog() 
Dim Count As Single 
Dim Dialog! As DialogSheet 
Dim Dialog? As DialogSheet 
Set Dialog! = DialogSheets("GetInitialData”) 
Set Dialog2 = DialogSheets("GetContactData") 


With Dialog] 

For Count = 1 To 8 
.EditBoxes(Count).Text = "Enter Value” 

Next Count 

End With 

Dialog1.Show 

InitializeContacts 

With Dialog2 
With .DropDowns(1) 


-RemoveAllltems 
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For Count = 1 To Parameters. NumberOfContacts 
_AddlItem Text:=Count, Index:=Count 
Next Count 
End With 
End With 
Dialog2.Show 
End Sub 


Sub SendInputToGraphsQ) 

Dim Count As Single 

Dim ProgSheet As Worksheet 

Set ProgSheet = Worksheets("Graphs") 

With ProgSheet 

For Count = 1 To 5 
.Cells(1, Count + 5) = "Contact " & Count 
.Cells(2, Count + 5) = Contacts(Count).Ct 
.Cells(3, Count + 5) = Contacts(Count).DMht 
.Cells(4, Count + 5) = Contacts(Count).Rh 
.Cells(5, Count + 5) = Contacts(Count).By 
.Cells(6, Count + 5) = Contacts(Count).CtSigma 
.Cells(7, Count + 5) = Contacts(Count).DMhtSigma 
.Cells(8, Count + 5) = Contacts(Count).RhSigma 
.Cells(9, Count + 5) = Contacts(Count).BySigma 
Next Count 

.Cells(2, 3) = Parameters. NumberOfContacts 

.Cells(3, 3) = Parameters. AcceptableMinRange 

.Cells(4, 3) = Parameters.SafetyRange 

.Cells(5, 3) = Parameters. TimeOnCourse 

.Cells(6, 3) = Parameters. DMho 

.Cells(7, 3) = Parameters. DirectionOfSeas 

.Cells(8, 3) = Parameters. NumberOfTnials 

.Cells(9, 3) = Parameters.StepSize 

End With 

End Sub 


35 


Sub ContactList_ Change() 
Dim Count As Single 
Dim CurrentDialog As DialogSheet 
Set CurrentDialog = Application. ActiveDialog 
With CurrentDialog 
For Count = 1 To 8 
.EditBoxes(Count).Text = "Enter Value" 
Next Count 
End With 
End Sub 


Sub NumberOfContacts_Change() 
Dim CurrentDialog As DialogSheet 
Set CurrentDialog = Application. ActiveDialog 
If IsNumeric(CurrentDialog.EditBoxes(1).Text) Then 
Parameters. NumberOfContacts = CurrentDialog.EditBoxes(1).Text 
End If 
End Sub 


Sub AcceptableMinRange Change() 
Dim CurrentDialog As DialogSheet 
Set CurrentDialog = Application. ActiveDialog 
If IsNumeric(CurrentDialog.EditBoxes(2).Text) Then 
Parameters. AcceptableMinRange = CurrentDialog.EditBoxes(2).Text 
End If 
End Sub 


Sub SafetyRange Change() 
Dim CurrentDialog As DialogSheet 
Set CurrentDialog = Application. ActiveDialog 
If IsNumeric(CurrentDialog.EditBoxes(3).Text) Then 
Parameters.SafetyRange = CurrentDialog.EditBoxes(3).Text 
End If 
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End Sub 


Sub TimeOnCourse_Change() 
Dim CurrentDialog As DialogSheet 
Set CurrentDialog = Application. ActiveDialog 
If IsNumeric(CurrentDialog.EditBoxes(4).Text) Then 
Parameters. TimeOnCourse = CurrentDialog.EditBoxes(4).Text 
End If 
End Sub 


Sub DMho_Change() 
Dim CurrentDialog As DialogSheet 
Set CurrentDialog = Application. ActiveDialog 
If IsNumeric(CurrentDialog.EditBoxes(5).Text) Then 
Parameters. DMho = CurrentDialog.EditBoxes(5).Text 
End If 
End Sub 


Sub DirectionOfSeas_Change() 
Dim CurrentDialog As DialogSheet 
Set CurrentDialog = Application. ActiveDialog 
If IsNumeric(CurrentDialog.EditBoxes(6).Text) Then 
Parameters.DirectionOfSeas = CurrentDialog.EditBoxes(6).Text 
End If 
End Sub 


Sub NumberOfTnals_Change() 
Dim CurrentDialog As DialogSheet 
Set CurrentDialog = Application. ActiveDialog 
If IsNumeric(CurrentDialog.EditBoxes(7).Text) Then 
Parameters. NumberOfTnials = CurrentDialog.EditBoxes(7).Text 
End If 
End Sub 
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Sub StepSize_Change() 
Dim CurrentDialog As DialogSheet 
Set CurrentDialog = Application. ActiveDialog 
If IsNumeric(CurrentDialog.EditBoxes(8).Text) Then 
Parameters.StepSize = CurrentDialog.EditBoxes(8).Text 
End If 
End Sub 


Sub CourseBox_Change() 
Dim CurrentDialog As DialogSheet 
Set CurrentDialog = Application.ActiveDialog 
If IsNumeric(CurrentDialog.EditBoxes(1).Text) Then 
Contacts(CurrentDialog.DropDowns. Value).Ct = CurrentDialog.EditBoxes(1).Text 
End If 
End Sub 


Sub SpeedBox_Change() 
Dim CurrentDialog As DialogSheet 
Set CurrentDialog = Application. ActiveDialog 
If IsNumeric(CurrentDialog.EditBoxes(2).Text) Then 
Contacts(CurrentDialog.DropDowns. Value). DMht = CurrentDialog.EditBoxes(2).Text 
End If 
End Sub 


Sub RangeBox_Change(Q) 
Dim CurrentDialog As DialogSheet 
Set CurrentDialog = Application. ActiveDialog 
If IsNumeric(CurrentDialog.EditBoxes(3).Text) Then 
Contacts(CurrentDialog. DropDowns. Value).Rh = CurrentDialog.EditBoxes(3).Text 
End If 
End Sub 


Sub BearingBox_Change() 
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Dim CurrentDialog As DialogSheet | 
Set CurrentDialog = Application. ActiveDialog 
If IsNumeric(CurrentDialog.EditBoxes(4).Text) Then 
Contacts(CurrentDialog. DropDowns. Value).By = CurrentDialog.EditBoxes(4).Text 
End If 
End Sub 


Sub CourseSigmaBox_Change() 
Dim CurrentDialog As DialogSheet 
Set CurrentDialog = Application. ActiveDialog 
If IsNumeric(CurrentDialog.EditBoxes(5).Text) Then 
Contacts(CurrentDialog.DropDowns. Value).CtSigma = CurrentDialog.EditBoxes(5).Text 
End If 
End Sub 


Sub SpeedSigmaBox_ChangeQ) 
Dim CurrentDialog As DialogSheet 
Set CurrentDialog = Application. ActiveDialog 
If IsNumeric(CurrentDialog.EditBoxes(6).Text) Then 
Contacts(CurrentDialog.DropDowns. Value). DMhtSigma = CurrentDialog.EditBoxes(6).Text 
End If 
End Sub 


Sub RangeSigmaBox_Change() 
Dim CurrentDialog As DialogSheet 
Set CurrentDialog = Application. ActiveDialog 
If IsNumeric(CurrentDialog.EditBoxes(7).Text) Then 
Contacts(CurrentDialog.DropDowns. Value).RhSigma = CurrentDialog.EditBoxes(7).Text 
End If 
End Sub 


Sub BeaningSigmaBox_Change() 


Dim CurrentDialog As DialogSheet 
Set CurrentDialog = Application. ActiveDialog 
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If IsNumenic(CurrentDialog.EditBoxes(5).Text) Then 


Contacts(CurrentDialog. DropDowns. Value).BySigma = CurrentDialog.EditBoxes(8).Text 
End If 


End Sub 
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APPENDIX C. DIALOGS 
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APPENDIX D. VISUAL BASIC SOURCE CODE FOR SIMULATION 


Option Explicit 


Function GetX(Rh As Single, By As Single) 
Application. Volatile 
GetX = Rh * Cos(By) 


End Function 


Function GetY(Rh As Single, By As Single) 
Application. Volatile 
GetY = Rh * Sin(By) 

End Function 


Sub GetBins() 
Dim Count As Single 
Dim Bin(0 To 20) As Single 
Parameters.Bins(0) = Application.NormSInv(0.01) 
For Count = 1 To 19 
Parameters.Bins(Count) = Application.NormSInv(0.05 * Count) 
Next Count 
Parameters.Bins(20) = Application. NormSInv(0.99) 
End Sub 


Sub ConvertCoordinates(Co As Single, Parameters As Parameter) 
Application. Volatile 
CCRider.Xo = GetX(Parameters. DMho, (Co * Application.PiQ) / 180)) 
CCRider. Yo = GetY(Parameters. DMho, (Co * Application.Pi() / 180)) 
End Sub 


Sub ConvertToRadians() 
Dim Count As Single 
Application. Volatile 
For Count = 1 To Parameters. NumberOfContacts 
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Contacts(Count).By = Contacts(Count).By * Application.Pi() / 180 
Contacts(Count).BySigma = Contacts(Count).BySigma * Application.Pi() / 180 
Contacts(Count).Ct = Contacts(Count).Ct * Application.Pi() / 180 
Contacts(Count).CtSigma = Contacts(Count).CtSigma * Application.Pi(Q) / 180 
Next Count 
End Sub 


Function GetError(Sigma As Single) As Single 
Dim P1 As Single 
Dim Mean As Single 
Application. Volatile 
Mean = 0 'for all error distributions 
Randomize 
Pl =Rnd 
If Sigma <> 0 Then 
GetError = Parameters.Bins(Int(21 * P1)) * Sigma / 2 
Else 
GetError = 0 
End If 
End Function 


Function GetTimeToCPA(Xo As Single, Yo As Single, Xc As Single, Yc As Single, Xp As 
Single, Yp As Single) As Single 
Dim X As Single 
Dim Y As Single 
Dim xDMhr As Single 
Dim yDMhr As Single 
Application. Volatile 
xDMhr = (Xc - Xo) * 33.75 
yDMhr = (Yc - Yo) * 33.75 
If Application. And(xDMhr, yDMhr) <> 0 Then 
GetTimeToCPA = -((xDMhr * Xp + yDMhr * Yp) / (kxDMhr * xDMhr + yDMhr * yDMhr)) 
Else 
GetTimeToCPA = 0 


44 


End If 
End Function 


Function GetMinSquaredRange(Xo As Single, Yo As Single, Xc As Single, Yc As Single, Xp As 
Single, Yp As Single, Time As Single) As Single 
Application. Volatile 
If Time <= 0 Then 
GetMinSquaredRange = GetNextSquaredRange(Xo, Yo, Xc, Yc, Xp, Yp, 0) 
Elself Time > Parameters. TimeOnCourse Then 
GetMinSquaredRange = GetNextSquaredRange(Xo, Yo, Xc, Yc, Xp, Yp, 
Parameters. TimeOnCourse) 
Else 
GetMinSquaredRange = GetNextSquaredRange(Xo, Yo, Xc, Yc, Xp, Yp, Time) 
End If 
End Function 


Function GetNextSquaredRange(Xo As Single, Yo As Single, Xc As Single, Yc As Single, Xp 
As Single, Yp As Single, Time As Single) As Single 
Dim X As Single 
Dim Y As Single 
Application. Volatile 
X = Xp + ((Xc - Xo) * 33.75) * Time 
Y = Yp + ((Yc - Yo) * 33.75) * Time 
GetNextSquaredRange = X * X + Y* Y 
End Function 


Function Simulate(Parameters As Parameter, Contacts() As Contact) As Variant 
Dim Course As Single 
Dim Trial As Single 
Dim Count As Single 
Dim TrialRh As Single 
Dim TnialBy As Single 
Dim TrialDMht As Single 
Dim TnialCt As Single 
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Dim TrialXc As Single 
Dim TrialYc As Single 
Dim TrialXp As Single 
Dim TrialYp As Single 
Dim Time As Single 
Dim n As Single 
Dim k As Single 
Dim AbsoluteMinSquaredRange As Single 
Dim MOE(1 To 2, 0 To 359) As Single 
Dim MinSquaredRange() As Single 
ReDim MinSquaredRange(1 To Parameters. NumberOfContacts) As Single 
Dim DBy As Single 
Application. Volatile 
ConvertToRadians 
GetBins 
For Course = 0 To 359 Step Parameters.StepSize 
i 
k=0 
ConvertCoordinates Course, Parameters 
For Trial = 1 To Parameters. NumberOfTnals 
For Count = 1 To Parameters. NumberOfContacts 
TrialRh = Contacts(Count).Rh + GetError(Contacts(Count).RhSigma) 
TnialBy = Contacts(Count).By + GetError(Contacts(Count).BySigma) 
TrialDMht = Contacts(Count). DMht + GetError(Contacts(Count).DMhtSigma) 
TrialCt = Contacts(Count).Ct + GetError(Contacts(Count).CtSigma) 
TnialXc = GetX(TnialDMht, TrialCt) 
Trial Yc = GetY(TnalDMht, TnialCt) 
TnialXp = GetX(TnialRh, TrialBy) 
TnialYp = GetY(TrialRh, TrialBy) 
Time = GetTimeToCPA(CCRider. Xo, CCRider. Yo, TrialXc, TnialYc, TrialXp, TrialYp) 
MinSquaredRange(Count) = GetMinSquaredRange(CCRider.Xo, CCRider. Yo, TrialXc, 
TnialYc, TnalXp, TrialYp, Time) 
Next Count 
AbsoluteMinSquaredRange = Application.Min(MinSquaredRange()) 
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If AbsoluteMinSquaredRange >= (Parameters.AcceptableMinRange 
Parameters. AcceptableMinRange) Then 
n=nt+] 
End If 
If AbsoluteMinSquaredRange >= (Parameters.SafetyRange * Parameters.SafetyRange) Then 
k=k+1 
End If 
Next Tnial 
MOE(1, Course) = n / Parameters. NumberOfTrials 
MOE(2, Course) = k / Parameters. NumberOfTrials 
Next Course 
Simulate = MOEQ 


End Function 


Sub SendDataToGraphs(MOEDataSet As Variant) 
Dim Count As Single 
Dim GraphSheet As Worksheet 
Set GraphSheet = Worksheets("Graphs") 
Application. Volatile 
For Count = 1 To 360 
With GraphSheet 
.Cells(Count + 1, 17) = MOEDataSet(1, (Count - 1)) 
.Cells(Count + 1, 18) = MOEDataSet(2, (Count - 1)) 
End With © 
Next Count 
End Sub 
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APPENDIX E. VISUAL BASIC SOURCE CODE FOR MAIN MODULE 


Option Explicit 
Dim Data As Variant 
Sub Main(Q) 
PrepareDialog 
SendInputToGraphs 
Data = Simulate(Parameters, Contacts) 
SendDataToGraphs (Data) 
End Sub 
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